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INTRODUCTION 


This  report  describes  the  work  done  on  Contract  No.  N00 1 23-76-C- 
1371  for  the  period  1  August  1977  to  31  October  1977. 

The  overall  program  is  divided  into  three  areas.  These  are  1) 
basic  scattering  code  development,  2)  reflector  antenna  code  development 
and  3)  basic  theoretical  studies  to  support  the  first  two  areas.  The 
following  sections  describe  the  research  accomplished  in  each  of  the 
three  areas  mentioned  above. 


II.  PROGRAM  SCOPE 

The  scope  of  the  work  under  Contract  No.  N001 23-76-C-l 371  is  to 
develop  the  necessary  theory,  algorithms  and  computer  codes  for  simu¬ 
lating  antennas  at  UHF  and  above  in  a  complex  ship  environment.  A 
revised  milestone  chart  for  the  three  year  program  is  shown  in  Table  I. 
The  milestone  chart  is  revised  from  that  shown  in  the  Annual  Report 
1  August  1976  to  31  July  1977.  These  revisions  were  agreed  upon  at  a 
project  review  meeting  at  NOSC  as  discussed  below. 


Project  Review  Meeti rig 

A  project  review  meeting  was  held  on  5  October  1977  at  NOSC  with 
C.H.  Walter,  R.C.  Rudduck,  and  W.D.  Burnside  representing  Ohio  State 
and  J.W.  Rockway,  J.  Logan,  and  S.T.  Li  representing  NOSC.  In  this 
meeting  the  current  status  of  the  computer  codes  was  described  in  detail. 
Computed  results  were  shown  to  illustrate  the  powerful  features  of  both 
the  Basic  Scattering  Code  and  the  Reflector  Antenna  Code.  Both  codes 
have  great  potential  for  solving  electromagnetic  interference  problems 
in  shipboard  antenna  systems.  Furthermore,  the  codes  can  be  used  to 
predict  the  performance  of  future  shipboard  antenna  designs. 

At  the  project  review  meeting  the  planned  deliverables  and  their 
projected  completion  dates  were  discussed.  It  was  agreed  to  make  some 
revisions  in  the  deliverables  which  will  not  significantly  change  the 
capability  of  the  codes  and  their  documentation. 

The  major  result  of  the  revisions  is  that  all  deliverable  items 
will  be  completed  by  the  end  of  the  third  year.  This  was  made  possible 
by  including  the  documentation  for  the  computer  codes  in  only  two  code 
manuals  instead  of  four  separate  manuals. 

At  the  request  of  Dr.  Rockway  several  command  features  will  be 
added  to  the  delivered  codes.  The  far  field  scattering  code  will  in¬ 
clude  an  option  for  range  dependence.  The  reflector  antenna  code  will 
include  the  capability  to  change  the  reflector  geometry  or  the  feed 
pattern  independently  so  that  data  inputs  can  be  substantially  reduced 
in  some  cases.  Furthermore,  specific  examples  of  Navy  reflector  antennas 
will  be  included  in  the  user's  manual.  The  planned  deliverables  and 
thei r  projected  completion  dates  are  shown  in  Table  II. 
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TABLE  I.  Revised  Milestone  Chart  (10/05/77) 


TABLE  II  -  PLANNED  DELIVERABLES 


BASIC  SCATTERING  CODE 

Projected  Completion 
Date 

October  31,  1977 
November  30,  1977 

October  31,  1978 
October  31.  1978 
March  31,  1979 
July  31,  1979 


DELIVERABLES 

Flat  Plates,  Box  and 
Cylinder  (Independent) 

Computer  Program  Package  (Delivered 

Oct.  5,  1977) 

User  Manual  (Typed  Rough  Draft) 

Coupled  Solution  For  Flat 
Plates,  Box  and  Cylinder* 

Computer  Program  Package 
User  Manual 

Rough  Draft  of  Code  Manual 
Code  Manual 


Will  supercede  previous  code 


REFLECTOR  ANTENNA  CODE  DELIVERABLES 


Completion  Date 

January  31,  1978 
February  28,  1978 


October  31,  1978 
October  31,  1978 


April  30,  1979 
July  31,  1979 
July  31,  1979 
July  31,  1979 

Will  supercede  previous  two  codes 


Far  Field  Code 

Computer  Program  Package 
User  Manual  (Typed  Rough  Draft) 

Near  Field  Code 

Computer  Program  Package 
User  Manual 

Code  with  Scattering 
from  Feed  and  Supports** 

Rough  Draft  of  Code  Manual 
Code  Manual 

Computer  Program  Package 
User  Manual 
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III.  BASIC  SCATTERING  CODE  DEVELOPMENT 


The  purpose  of  this  section  is  to  describe  the  present  status  of 
the  basic  scattering  code  development  program  for  the  analysis  of  antennas 
in  a  complex  shipboard  environment.  The  Geometrical  Theory  of  Diffraction 
(GTD)  is  being  used  to  develop  algorithms  to  solve  for  the  scattering 
from  basic  plate  and  cylinder  structures.  These  simple  components  can 
then  be  combined  to  form  box-like  structures  with  a  nearby  finite  ellip¬ 
tic  cylinder  that  can  represent  the  various  component  structures  of  a 
ship.  The  algorithms  are  being  implemented  into  a  user-oriented  computer 
code. 


In  this  period,  the  NEC  code  has  been  interfaced  with  the  GTD  scat¬ 
tering  code.  The  GTD  code  will  accept  as  input  data  the  output  data  of 
the  NEC  code  which  describes  wire  segment  location,  orientation,  and 
excitation.  The  GTD  code  is  written  at  present  to  use  only  the  wire 
segment  capabilities  of  the  NEC  code  and  not  the  patch  segments.  In 
order  to  simplify  the  interfacing  of  the  two  codes  and  to  make  the  GTD 
code  as  fast  and  efficient  as  possible,  the  NEC  currents  are  assumed  to 
be  the  value  of  the  current  at  the  center  of  a  segment  with  a  cosine 
current  distribution.  This,  of  course,  is  not  the  same  distribution  as 
used  in  the  NEC  code  but  as  is  well  known  the  radiation  pattern  is 
insensitive  to  the  current  distribution  for  segments  shorter  than  >/4. 

The  GTD  code  is  normalized  in  the  far  field  to  exp(-jkR)/R.  The  fields 
are  output  in  volts/meter,  however,  it  should  be  noted  that  unless  a 
range  R  is  defined  the  units  are  actually  in  volts.  The  directive  gain 
is  also  calculated  by  using  the  power  radiated  as  supplied  by  the  NEC 
code. 


The  interfacing  of  the  NEC  and  GTD  codes  has  been  tested  by  using 
two  examples  obtained  from  Mr.  Jim  Logan  at  NOSC.  Both  examples  are 
composed  of  antennas  that  were  modeled  by  wire  segments  over  a  perfectly 
conducting  infinite  ground  plane.  Using  the  assumption  that  the  current 
distributions  obtained  in  this  way  will  not  be  perturbed  appreciably  by 
slightly  changing  the  environment,  the  currents  from  the  NEC  code  can 
be  used  to  simulate  the  radiation  pattern  for  the  antennas  mounted  over 
a  finite  plate  which  is  over  an  infinite  ground  plane  as  shown  in 
Figure  1 . 

To  test  the  compatibility  of  the  current  distribution  used  in  the 
GTD  code  with  the  current  weights  obtained  from  the  NEC  code,  the  GTD 
cede  was  used  to  calculate  the  field  for  the  example  problem  of  four 
quarter-wavelength  dipoles.  The  locations  of  the  four  dipoles  are 
illustrated  in  Figure  1  as  the  four  solid  lines  parallel  to  the  x-axis. 
The  infinite  ground  plane  is  coincident  with  the  square  plate  for  this 
first  test.  The  dipoles  were  divided  into  a  total  of  20  segments.  The 
fields  obtained  from  this  calculation  should  have  exactly  duplicated 
the  fields  as  calculated  by  the  NEC  code.  However,  it  was  found  that 
the  GTD  code  was  low  by  a  factor  of  approximately  1.55.  This  is  very 
close  to  a  factor  of  one-half  pi  which  indicated  the  GTD  code's  cosine 
current  distribution  does  not  give  the  correct  level  when  used  for  a 


small  segment  size.  The  pattern  shapes  for  the  two  results,  however, 
are  identical.  When  a  factor  of  one-half  pi  was  added  to  the  GTD  results 
the  directive  gains  were  within  0.2  dB  as  shown  in  Figures  2a  and  2b  for 
the  plane  4>=0°  and  $=90°,  respectively.  It  should  be  noted  that  when 
the  four  quarter-wavel ength  dipoles  with  a  cosine  distribution  weighted 
to  the  center  current  value  is  used  directly  in  the  GTD  code,  the  field 
is  very  close  to  the  NEC  result  and  therefore  does  not  need  the  factor 
of  one-half  pi.  The  current  distribution  used  in  the  GTD  code,  there¬ 
fore,  is  only  off  for  small  sized  segments.  The  factor  of  one-half  pi 
has  been  included  in  the  code  only  when  the  NEC  option  is  used  and 
should  be  valid  if  the  segment  lengths  used  are  on  the  order  of  approxi¬ 
mately  0.1  wavelength  or  less. 

The  GTD  code  was  next  used  to  show  how  the  directive  gain  changes 
when  the  four  dipoles  are  placed  over  a  square  plate  over  a  ground  plane 
as  shown  in  Figure  1.  The  comparisons  between  the  infinite  ground  plane 
and  square  plate  cases  are  shown  in  Figures  3a  and  3b  for  the  two  pattern 
planes. 

In  order  to  confirm  the  above  conclusions  on  the  current  distribu¬ 
tions,  a  second  example  was  used.  This  one  is  composed  of  four  sets  of 
crossed  dipoles  as  shown  in  Figure  1  where  solid  and  a  dashed  lines  are 
used  to  represent  the  crossed  dipoles.  The  results  obtained  from  NOSC 
for  this  example  were  obtained  from  an  AMP  code  print  out  rather  than  a 
NEC  code  print  out;  however,  the  results  are  still  useful  for  comparison 
purposes.  The  GTD  and  AMP  results  for  48  segment  dipoles  over  an  infin¬ 
ite  ground  plane  are  compared  in  Figures  4a  and  4b,  where  Figure  4a  shows 
the  vertical  component  and  Figure  4b  the  horizontal  component.  The 
solid  curves  represent  the  GTD  code  with  a  factor  of  one-half  pi  added 
and  the  X's  show  the  AMP  result.  Again  the  two  results  compare  within 
0.2  dB.  Square  plate  results  are  also  shown  in  Figures  4a  and  4b. 
Directive  gains  for  the  infinite  ground  plane  and  square  plate  cases 
are  compared  in  Figures  5a  and  5b  for  the  4>=0°  and  <j>=45°  planes, 
respectively. 

The  second  example  confirms  the  conclusion  of  the  first  example 
that  the  GTD  code  and  NEC  code  appear  now  to  be  properly  interfaced. 

It  should  be  noted  that  the  directive  gain  is  being  computed  using  the 
power  radiated  as  given  in  the  NEC  code.  If  this  value  is  off  in  the 
NEC  code  then  the  GTD  directive  gain  will  be  off  by  the  same  factor 
with  everything  else  being  the  same.  An  alternative  method  of  computing 
the  directive  gain  would  be  to  integrate  the  volumetric  pattern  of  the 
GTD  code  and  use  that  as  the  power  radiated.  This  is  possible  but  would 
not  be  an  efficient  method  since  it  would  be  very  costly  to  compute  all 
the  needed  pattern  values.  This  is  especially  true  for  the  NEC  code 
interface,  since  in  general  a  large  number  of  current  segments  will  be 
used  which  increases  the  computer  time  needed  for  each  pattern. 

The  basic  scattering  code  that  allows  simulation  of  concave  and 
convex  plate  structures  and  a  finite  elliptic  cylinder  structure  with¬ 
out  coupling  between  the  plates  and  cylinder  was  completed  and  delivered 
to  NOSC  on  5  October  1977.  The  NEC  code  was  interfaced  in  this  code. 


Comparison  of  the  directive  gain  of  four  dipoles  over  an 
infinite  ground  plane  obtained  from  the  NEC  '■ode  and  from 
the  GTD  code  using  the  NEC  currents  (<J>=0°  vertical 
pol arization) . 


DIRECTIVE  GAIN 


DIRECTIVE  GAIN  (dB) 


Figure  3b.  Comparison  or  the  directive  gain  of  four  dipoles  over  an 
infinite  oround  with  four  dipoles  over  a  squate  plate 
over  an  in,  mite  ground  plane  (*>=90°,  horizontal 
polarization) . 


DIRECTIVE  GAIN 


DIRECTIVE  GAIN 


Auor  ,NpINITE  GROUND 

^  *  M  M  PI 

_  SQUARE  PLATE  OVER 
INFINITE  GROUND 


Figure  5a.  Comparison  of  the  total  directive  gain  of  four  sets  of 

crossed  dipoles  over  an  infinite  ground  and  over  a  square 
plate  over  an  infinite  ground  U=0°). 
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Figure  5b.  Comparison  of  the  total  directive  gain  of  four  sets  of 

crossed  dipoles  over  an  infinite  ground  and  over  a  square 
plate  over  an  infinite  ground  (<t>=45°). 


but  the  one-half  pi  factor  and  improvements  in  the  print  out  representa¬ 
tion  of  the  directive  gain  have  been  added  since  delivery.  An  option 
allowing  a  range  in  the  far  field  to  be  specified  has  also  been  added. 
These  updates  can  easily  be  implemented  in  the  delivered  code.  The  user's 
manual  is  also  under  preparation  and  should  be  completed  early  in  December 
1977.  The  user's  manual  will  explain  in  detail  the  input  parameters  and 
output  data  of  the  code.  Specific  examples  will  illustrate  the  versatil¬ 
ity  of  the  code  and  will  clarify  many  of  the  subtleties  involved  in  using 
it  to  simulate  scattering  structures. 


IV.  REFLECTOR  ANTENNA  CODE  DEVELOPMENT 

The  major  task  of  the  present  effort  under  the  reflector  antenna 
code  development  is  to  develop  a  user-oriented  computer  program  package 
by  which  the  far  field  pattern  of  a  typical  Navy  reflector  antenna  can 
be  calculated.  Feed  blockage  and  scattering  effects  are  not  included 
in  this  phase.  These  will  be  included  later  (see  Table  I). 

The  theoretical  approach  for  computing  the  far  field  pattern  of 
the  general  reflector  is  based  on  a  combination  of  the  Geometrical  Theory 
of  Diffraction  (GTD)  and  Aperture  Integration  (AI)  techniques.  AI  will 
be  used  to  compute  the  main  beam  and  near  sidelobes;  GTD  will  be  used  to 
compute  the  wide-angle  sidelobes  and  the  backlobes.  To  implement  the 
computer  algorithms  based  on  these  theories,  efficient  ways  are  being 
developed  to  handle  calculations  involving  the  feed  pattern  and  the 
reflector  geometry.  The  geometry  of  the  reflector  rim  is  treated  as 
piece-wise  linear. 

Since  the  majority  of  Navy  reflector  antennas  have  parabolic  sur¬ 
faces,  only  the  class  of  parabolic  surfaces  will  be  implemented  in  the 
computer  code.  The  code  for  the  reflector  geometry  will  be  flexible 
enough  to  include  off-set  fed  reflectors  and  general  reflector  rim 
shapes  such  as  elliptical  and  rectangular  with  chopped  corners. 

An  efficient  way  to  handle  AI  is  necessary  because  of  its  time- 
consuming  nature  for  electrically  large  reflectors.  Fortunately,  AI  is 
needed  only  for  the  main  beam  and  near  sidelobe  regions  of  the  pattern. 
GTD  is  very  efficient  for  calculating  most  of  the  pattern.  Two  techni¬ 
ques  are  being  implemented  to  provide  efficient  AI:  the  polynomial -fi t 
method  for  the  feed  pattern  and  the  rotating  grid  method  for  the  numeri¬ 
cal  integration  required  to  calculate  the  far  field  pattern. 


A.  Feed  Pattern 

The  polynomial -fi t  method  is  used  to  obtain  an  analytic  fit  for 
each  measured  pattern  cut  of  the  volumetric  feed  pattern.  This  method 
provides  a  computationally  efficient  way  of  calculating  the  aperture 
field  without  requiring  large  amounts  of  computer  storage  for  the  mea¬ 
sured  feed  pattern.  Only  relatively  few  coefficients  need  to  be  stored 
for  essentially  complete  feed  pattern  information.  Furthermore,  the 


polynomial-fit  method  has  the  advantages  of  flexibility  and  simplicity 
for  general  feed  patterns.  No  cut-and-try  procedures  are  needed;  the 
polynomial  coefficients  can  be  computed  automatically  from  the  measured 
feed  pattern  input. 

Further  improvement  has  been  made  for  fitting  measured  feed  patterns 
in  the  spillover  region,  i.e.,  where  the  feed  radiation  is  not  reflected 
by  the  reflector.  At  the  low  field  levels  encountered  in  the  spillover 
region,  a  piecewise  linear  approximation  was  found  to  be  a  better  repre¬ 
sentation  for  this  part  of  the  feed  pattern.  Thus  the  computer  sub¬ 
routines  were  modified  to  use  a  combination  of  polynomial -fi t  for  the 
major  part  of  the  feed  pattern  and  a  piecewise  linear  fit  for  most  of 
the  spillover  region  of  the  pattern. 

Several  subroutines  have  been  developed  to  treat  a  feed  with  a 
general  pattern  and  general  polarization.  These  subroutines  decompose 
the  feed  radiation  into  components  to  provide  the  necessary  information 
for  the  following:  1)  aperture  fields,  2)  GTD  edge  illuminations,  and 
3)  direct  feed  radiation. 

In  the  computer  code  the  measured  feed  pattern  is  inputted  at  the 
angles  <}>=<t>n>  n=l,2, — N^.  The  theory  for  obtaining  the  combined  poly¬ 
nomial-linear  fit  for  a  specific  input  feed  pattern  cut  at  <t>n  is  given 
in  the  Annual  Report  [1].  The  feed  pattern  in  planes  other  than  those 
corresponding  to  the  input  pattern  cuts  is  calculated  by  linear  inter¬ 
polation.  As  a  first  step  the  nearest  input  pattern  angles  4>q  and  $p 
are  determined  for  the  desired  feed  pattern  angle  <t>  as  shown  in  Figure  6. 
Then  the  pattern  is  calculated  from 


=  Qdp  +  O-Qop)  0) 

where 

%P  =  4^  (2) 


Various  states  of  feed  polarizations  can  be  treated  by  using  two 
orthogonal  linear  polarizations.  The  two  linear  polarization  components 
are  treated  as  those  corresponding  to  dipoles  oriented  along  the  x-axis 
and  the  y-axis,  respectively.  Thus  the  electric  field  of  the  feed  is 
given  by 


E  [-t’x  AX  fx('^»'l>  )  “  Ay  )] 


y  y  yv 


-jkR' 

R' 


(3) 


where  -ox  and  -Gy  are  the  unit  vectors  corresponding  to  the  respective 
dipole  polarizations.  The  modified  spherical  coordinate  system  ( R ' , ^ ' ) 


Figure  6.  Interpolation  of  the  feed  pattern 
in  the  41-plane. 


for  the  feed  pattern  is  shown  in  Figure  7a.  The  unit  vector  ey  is  defined 
by  the  polar  angle  0y  which  js  measured  from  the  y-axis  as  shown  in 
Figure  7b.  The  unit  vector  0X  is  defined  by  the  angle  0X  which  is  mea¬ 
sured  in  the  same  way  but  from  the  x-axis. 


y 


y 


Figure  7.  Coordinate  systems  for  the  primary  feed. 


18 


Each  pattern  component  fx  and  fy  must  be  consistent  for  the  direction 
*=0;  thus  fx (0»4» ' )  and  f y ( 0 » 4> 1 ;  must  be  constant  for  all  values  *'.  The 
feed  polarization  is  controlled  by  the  complex  constants  Ax  and  A  For 

example,  circular  polarization  can  be  obtained  by  setting  Ay=ijAx. 

The  unit  vector  ex  can  be  expressed  in  terms  of  e  and  *  components 
of  the  standard  spherical  coordinate  system  (R,e,*): 


e  =  I—  [-0  cose  cos*  + 
Bx 

=  5—  [e  cos*  cos*'  + 

by 


where  *=tt-0  and  *'=*  are  the 


*  sin*] 

*  sin*']  (4) 

coordinate  angles  for  the  feed  pattern  and 


Bx  =  J l-sin^*  cos^*‘ 

Similarly  the  unit  vector  0y  is  given  by 

A  "j  A  A 

0  =  ~  [e  cos*  sin*'  -  *  cos*'] 

*  ^\/ 


(5) 


(6) 


where 


B  =  Jl-sin1-*  sin 


2.,  „,-„2 , 1 


(7) 


Thus  the  feed  pattern  can  be  expressed  as 

.  -JkR' 

FT 


E  ■  -  [e  f„  +  ♦  f+]  -  w 


(8) 


where 


Ax  A 

f0  =  cos*  cos*'  f  (*,<!»')  +  cos*  sin*'  ^  f  (*,*') 
x  y 


(9) 


and 


A  A 

sin*'  )  -  cos*'  n^f  (*,*’) 

Dx  x  y  y 


(10) 
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Equation  (8)  is  used  to  calculate  the  feed  pattern  for  the  direct 
feed  radiation  and  also  to  determine  the  aperture  field.  For  a  para¬ 
boloidal  surface  with  a  feed  at  the  focus  the  aperture  field  at  point 
P(x,y,z)  as  shown  in  Figure  8  is  given  by 


e-JkRo 


Ex  =  -  [cos4> 1  fgU.*')  +  sin<j>'  f^U.^'  )F-^r 


and 

a  e"jkRo 

E  =  -  [sin*1  f0(«|»,«t> 1 )  -  cos<j> 1  f^ <t> ' )] — pi — 


(12) 


The  necessary  coordinate  information  for  the  aperture  field  is  defined 
by  the  point  of  reflection  on  the  reflector  surface  with  coordinates 
x,y  as  shown  in  Figure  8.  Thus 


(14) 


<j/  =  tan“l  £_  (15) 

ZP 

and 

<p'  =  tan"^  X-  •  (16) 

x 

Rq  is  the  pathlength  from  the  focus  to  the  aperture  plane  which  is  defined 
as  the  plane  perpendicular  to  the  z-axis  and  passing  through  the  reference 
point  P0(x0,y0,z0)  on  the  reflector  rim.  Thus 
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The  field  of  the  feed  incident  on  the  reflector  rim  is  needed  to 
calculate  the  diffracted  field  using  the  GTD  approach.  For  this  purpose 
it  is  convenient  to  decompose  the  feed  polarization  vectors  ex  and  9 
into  rectangular  components.  Thus  y 


9x  =  ~x  Bv  +  y 


sin^  cos*1  sin*1  '  sin*  cos*  cos*1 


2 

e  =  x  y  cos  * 1  sin<p 1  _  “  B  “  sin*  cost  sin*1  .  . 

-y  B  y  y  Bw  ^  ' 

7  y  7  y 

where  B„  and  By  are  given  by  Equations  (5)  and  (7).  Thus  the  field  of 
the  feed  is  expressed  in  rectangular  components  by 

f  A  1  - jkR 1 

Ev  =  "Bx  Ax  fx^’^'^  +  sin^  cos*'  sin$'  U'M')  —  (22) 


B  y' 

y  7 


i  2  Ax  e'jkR' 

Ey  =  sin  *  c o s 4> 1  s i n 4> 1  ^  f x ( ^ , <*> ’ )  -  By  Ay  fy(*,*')  —  (23) 

X 


r  a  A  -jkR1 

=  -sin*  cos*  cos* 1  jp  fx(*,*' )  +  sin*'  fy(*,*' )  —  .  (24 


Rotating  Grid  Method 


The  rotating  grid  provides  greatly  improved  efficiency  for  the  far 
field  pattern  computations.  With  this  approach  the  y-integration  is 
carried  out  for  each  column  of  the  aperture  and  each  one-dimensional 
integration  result  is  stored.  The  stored  values  for  the  y-integration 
are  then  used  for  each  pattern  angle  e  in  the  plane  perpendicular  to 
the  y-axis;  thus  the  efficiency  approaches  that  of  a  one-dimensional 
integration.  Even  though  the  integration  grid  must  be  rotated  to  obtain 
the  pattern  in  other  planes,  the  required  grid  rotation  is  computation¬ 
ally  much  faster  than  the  numerous  two-dimensional  integrations  that 
would  otherwise  be  required.  The  computer  algorithms  for  the  rotating 
grid  method  have  been  developed  and  the  method  has  been  programmed  on 
the  OSU  Datacraft  computer. 


Some  typical  far  field  pattern  results  from  the  rotating  grid 
method  were  shown  in  the  Annual  Report  [1].  The  accuracy  of  these  results 
were  checked  by  using  the  rectangular  subaperture  method  described  in  the 
Third  Quarterly  Report  [2].  It  was  found  that  very  accurate  results 
could  be  achieved  with  very  good  computational  efficiency.  The  rotating 
grid  method  was  found  to  run  7-10  times  faster  than  the  rectangular  sub¬ 
aperture  method  as  shown  in  Table  III. 


C.  GTD  Method 

The  GTD  will  be  used  to  calculate  the  wide-angle  sidelobe  and  back- 
lobe  regions  for  the  reflector  antenna  computer  code.  Work  is  continuing 
on  the  development  of  the  GTD  portion  of  this  code.  Advantage  is  being 
taken  of  the  existing  GTD  computer  subroutines  for  the  edge  diffraction 
functions  and  the  corner  diffraction  functions. 

Part  of  the  geometry  subroutine  from  the  flat  plate  scattering  code 
has  been  modified  for  use  in  the  GTD  portion  of  the  reflector  antenna 
code.  The  3-D  geometry  of  the  reflector  rim  is  determined  from  the  inter¬ 
section  of  the  projected  rim  shape  and  the  parabolic  reflector  surface. 

The  projected  rim  shape  is  specified  as  piece-wise  linear  as  shown  in 
Figure  9.  The  input  for  the  rim  geometry  is  specified  by  the  X  and  Y 

coordinates  of  each  junction  point  on  the  rim.  The  Z  coordinate  of  each 

junction  is  given  by 

X2  +  V2  -  D2 

Z  »  - te - 2  (25) 


where  F  is  the  focal  length  of  the  parabola  and  pQ  is  the  radius  to  the 
reference  point  P  which  defines  the  aperture  plane  and  is  given  by 
Equation  (18). 

The  geometry  associated  with  the  reflector  rim  is  calculated  and 
stored  for  use  in  calculating  all  of  the  diffracted  fields  via  the  GTD. 
The  following  information  for  the  geometry  is  stored  for  each  linear 
segment  of  the  reflector  rim: 

1)  the  edge  unit  vectors, 

2)  the  unit  vector  normal  to  the  parabolic  surface, 

3)  the  unit  binormal  vector  (perpendicular  to  the  edge 
and  normal  vectors),  and 

4)  the  bounds  on  the  permissible  range  for  the  diffraction 
angle. 
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Figure  9.  Reflector  rim  geometry. 


In  order  to  improve  the  computational  efficiency  of  the  code  a 
shadow  boundary  subroutine  has  been  developed  to  calculate  the  angles 
at  which  the  shadow  boundaries  occur  in  the  plane  in  which  the  pattern 
is  calculated.  Thus  the  subroutine  will  determine  the  pattern  angles 
for  which  the  direct  feed  pattern  (feed  spillover)  must  be  included  be¬ 
fore  the  pattern  is  computed.  This  will  eliminate  the  need  to  check 
whether  the  source  or  feed  must  be  included  for  each  pattern  angle,  as 
must  be  done  in  scattering  codes  with  more  than  one  scatterer. 

When  the  GTD  method  is  used  the  reflector  code  will  increment 
around  the  rim  and  determine  whether  a  diffraction  occurs  for  each  linear 
rim  segment.  This  will  be  done  by  comparing  the  diffraction  angle  with 
the  bounds  on  the  permissible  range  of  angles.  If  no  diffraction  occurs 
for  that  segment  the  code  will  check  the  next  rim  segment.  If  a  dif¬ 
fraction  does  occur,  the  diffraction  point  and  the  vector  for  the  inci¬ 
dent  ray  from  the  feed  will  be  calculated.  This  procedure  is  the  same 
as  that  used  for  the  flat  plate  scattering  code  except  that  the  geometry 
information  associated  with  the  parabolic  reflector  surface  is  changed. 

The  incident  field  at  the  diffraction  point  is  calculated  from 
the  feed  pattern  subroutines  discussed  in  a  previous  section.  The 
incident  field  is  decomposed  into  rectangular  components  by  use  of 
Equations  (22)  through  (24).  The  code  then  calculates  the  components 


of  the  incident  field  perpendicular  and  parallel  to  the  rim  segment. 

The  edge  diffraction  coefficients  and  the  corner  diffraction  coeffici¬ 
ents  for  the  respective  perpendicular  and  parallel  components  are  cal¬ 
culated  in  the  same  way  as  for  the  flat  plate.  See  Reference  [3]. 
Finally  the  diffracted  fields  are  calculated  and  summed  for  each  rim 
segment. 

Results  have  been  obtained  from  the  preliminary  GTD  code  and  it 
was  found  that  the  reflector  rim  must  be  subdivided  into  many  straight 
line  segments  because  of  its  3-D  curvature.  This  is  true  even  though 
the  projected  rim  shape  may  be  composed  of  only  a  few  straight  line  seg 
ments,  e.g.,  the  rectangular  aperture.  A  criterion  for  subdividing  the 
reflector  rim  is  currently  being  checked. 


V.  THEORETICAL  STUDIES 

The  completion  date  for  the  slope  diffraction  report  has  been 
changed  to  31  January  1978.  This  delay  has  been  caused  by  a  shift  of 
emphasis  to  the  work  on  double  diffraction  which  was  being  carried  out 
by  Professor  Tiberio,  a  visiting  professor  from  the  University  of 
Florence,  Italy.  Professor  Tiberio  has  been  forced  to  return  to  Italy 
sooner  than  expected  and  Professor  Kouyoumjian  had  to  drop  the  slope 
diffraction  work  in  order  to  help  Professor  Tiberio  complete  his  phase 
of  the  double  diffraction  work.  It  is  anticipated  that  the  preparation 
of  the  report  on  slope  diffraction  will  be  resumed  in  November  or  by 
December  at  the  latest  and  completed  by  31  January  1978. 
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